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Abstract  -  Harbor  sun’eillance  above  and  below  the  sea 
surface  depends  on  sensors  such  as  sun’eillance  radar  and 
multibeam  sonar.  These  sensors  attempt  to  detect  and  track 
moderately  obsen’able  targets  such  as  small  boats  or  hu¬ 
man  divers  in  environments  which  often  are  characterized  by 
heavy-tailed  backgrounds.  Target  tracking  in  heavy-tailed 
environments  is  challenging  even  for  moderately  strong  tar¬ 
gets  due  to  the  more  frequent  occurrences  of  target-like  out¬ 
liers.  One  strategy  for  increased  robustness  is  to  use  the 
backscattered  signal  strengths  together  with  the  kinematic 
measurements  in  the  tracking  method.  This  paper  proposes 
two  new  amplitude  likelihoods  for  target  tracking  in  heavy- 
tailed  backgrounds.  The  first  likelihood  works  by  incor¬ 
porating  the  uncertainty  of  the  background  estimate.  The 
second  likelihood  explicitly  treats  the  background  as  heavy¬ 
tailed  using  the  K -distribution. 

Keywords:  Data  association,  Detection,  Feature  Aided 
Tracking 

1  Introduction 

In  harbor  surveillance  sensors  such  as  radars  and  sonars  are 
used  to  detect  and  track  various  targets.  The  sensor  pro¬ 
vides  an  array  of  observations  to  be  further  analyzed  through 
signal  processing  and  data  analysis.  After  preprocessing, 
such  as  matched  filtering,  a  detector  is  employed  to  extract 
point  measurements  from  the  array.  This  is  done  by  setting  a 
threshold  according  to  a  Constant  False  Alarm  Rate  (CFAR) 
criterium  [9],  Then  tracks  are  established  and  maintained 
by  feeding  these  exttacted  point  measurements  to  a  target 
tracking  method. 

The  targets  we  have  in  mind  are  small  boats  in  radar  data 
or  human  divers  in  sonar  data.  Both  can  be  expected  to  have 
a  Signal  to  Noise  Ratio  (SNR)  somewhere  between  OdB  and 
20dB.  This  means  that  if  the  target  is  to  be  detected  with 
a  reasonably  high  probability  one  must  allow  some  false 
alarms  to  occur.  In  order  to  establish  and  maintain  a  track 
on  such  a  target  one  must  determine  which  detections  are 
actually  from  the  target,  and  which  detections  are  to  be  con¬ 
sidered  as  false  alarms  or  clutter. 


A  popular  solution  to  this  problem  of  data  association  is 
found  in  the  Probabilistic  Data  Association  Filter  (PDAF) 
and  its  multi-target  version  Joint  Probabilistic  Data  Asso¬ 
ciation  Filter  (JPDAF)  [4],  One  way  to  improve  the  per¬ 
formance  of  the  PDAF  is  the  utilization  of  Amplitude  In¬ 
formation  (AI).  In  [11]  the  PDAFAI  (PDAF  with  Ampli¬ 
tude  Information)  was  therefore  presented.  While  the  PDAF 
only  uses  kinematic  measurements  zk(i ),  the  PDAFAI  also 
uses  the  corresponding  amplitudes  ak  (i)  in  order  to  decide 
which  measurement  is  most  likely  to  originate  from  the  tar¬ 
get.  Of  course  the  utilization  of  AI  can  be  incorporated  in 
other  tracking  methods  as  well.  Although  Al-based  tracking 
methods  have  been  around  for  some  time,  the  usage  of  AI 
does  not  appear  to  be  commonplace  in  practical  systems. 

Figure  1  illustrates  the  difference  between  the  conven¬ 
tional  PDAF  and  the  PDAFAI.  It  is  meant  to  highlight  the 
fact  that  the  two  methods  primarily  differ  in  the  data  associ¬ 
ation  part,  and  not  in  the  state  estimation  part.  Both  meth¬ 
ods  calculate  a  conditional  state  estimate  xk\ k{i)  for  each 
measurement  based  on  purely  kinematic  information.  The 
PDAF  solely  depends  on  this  kinematic  information  in  the 
process  of  lumping  these  conditional  state  estimates  together 
to  the  final  a  posteriori  state  estimate  xk\k.  The  PDAFAI 
uses  the  amplitudes  of  the  measurements  as  well  in  this  pro¬ 
cess.  Since  the  PDAFAI  has  access  to  more  information  than 
the  PDAF,  it  can  be  expected  to  outperform  the  PDAF,  as 
long  as  the  background  and  target  amplitude  statistics  are 
adequately  modeled. 

Most  work  on  AI  assumes  that  the  background  is  Gaus¬ 
sian  or  Rayleigh  distributed  [11,14]  with  a  known  power. 
Such  assumptions  may  or  may  not  be  adequate  in  the  real 
world.  Experimental  evidence  [1]  has  indicated  that  more 
heavy-tailed  background  models  should  be  considered.  In 
particular,  the  /^-distribution  [15]  is  popular  in  the  radar  and 
sonar  communities. 

Many  papers  have  been  published  on  detectors  in  heavy¬ 
tailed  clutter,  see  for  example  [3]  or  [8].  However,  such  pa¬ 
pers  make  no  mention  of  target  tracking.  This  is  problem¬ 
atic,  as  a  detection  system  in  heavy-tailed  noise  may  need  to 
operate  at  high  false  alarm  rates  («  10-2).  The  decision  as 
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Figure  1:  Schematic  description  of  the  PDAF  and  PDAFAI  methods.  The  dashed  connection  represents  the  flow  of  ampli¬ 
tude  information  from  detector  to  tracker,  which  only  is  present  in  the  PDAFAI. 


to  whether  the  target  is  present  or  not  can  obviously  not  be 
made  based  on  the  detection  process  alone. 

One  way  to  overcome  this  challenge  is  to  integrate  the 
data  before  the  detection  decision  is  made,  as  done  in  so- 
called  Track-Before-Detect  (TBD)  [5,14],  For  very  low 
SNR,  when  the  target  return  cannot  be  expected  to  exceed 
any  reasonable  threshold,  this  is  the  only  feasible  approach 
to  target  tracking.  Otherwise  TBD  is  not  worthwhile  the 
very  large  computational  expenses. 

In  [5]  it  was  observed  that  effects  similar  to  the  well 
known  phenomena  of  CFAR  loss  [9]  could  deteriorate  the 
performance  of  TBD  methods.  Since  TBD  only  relies  on 
amplitude  information,  it  is  more  sensitive  to  the  interpre¬ 
tation  of  the  background  than  conventional  tracking  meth¬ 
ods.  To  prevent  the  method  from  being  misled  by  misin¬ 
terpretation  of  the  background  it  was  necessary  to  make  its 
amplitude  likelihood  more  conservative.  This  was  done  by 
marginalizing  over  the  true,  but  unknown,  noise  parameter. 

This  paper  presents  further  work  on  conservative  ampli¬ 
tude  likelihoods.  While  [5]  focused  on  TBD  methods,  this 
paper  focuses  on  more  traditional  PDAF  and  PDAFAI  meth¬ 
ods.  Analogous  to  what  was  done  in  [5]  we  first  develop 
a  proper  treatment  of  uncertain,  but  Rayleigh  distributed, 
background  noise.  Then  we  establish  the  treatment  of  am¬ 
plitude  information  in  the  more  challenging  case  of  heavy¬ 
tailed  clutter. 

The  paper  is  organized  as  follows:  Section  2  presents  the 
framework  in  which  our  tracking  methods  are  developed.  In 
Section  3  we  describe  the  tracking  methods,  with  particu¬ 
lar  attention  to  the  amplitude  likelihoods.  We  also  describe 
the  detection  process  preceding  target  tracking,  and  how  the 
various  quantities  involved  are  estimated.  A  conclusion  is 
given  in  Section  4.  Details  regarding  the  numerical  evalua¬ 
tion  of  the  amplitude  likelihoods  are  left  for  the  Appendix. 
Simulation  results  are  published  in  a  separate  paper  [6], 

2  Conceptual  framework 

In  this  paper  we  restrict  attention  to  the  single-target  track¬ 
ing  problem.  Although  the  real  world  most  often  poses 
problems  that  are  suitably  viewed  as  multiple  target  prob¬ 
lems,  the  single-target  problem  must  be  properly  understood 
before  one  can  move  onto  multiple  target  problems.  The 
utilization  of  AI  is  one  such  challenge  where  significant 
progress  can  and  should  be  made  in  the  single-target  case. 


The  single-target  tracking  problem  can  be  phrased  in 
terms  of  the  optimal  Bayes  equations, 

p{xk\Zk~l)  =  J  p(xk\xk-i)p(xk-1\Zk~1)dxk-i, 

p(xk\Zk)  oc  p{Zk\xk)p{xk\Zk~1).  (1) 

Our  objective  is  to  evaluate  the  posterior  Probability  Density 
Function  (PDF)  p(xk  \  Zk)  from  the  set  of  received  measure¬ 
ments  Zk  =  {Zi, . . . ,  Zkj  where  each  Zk  contains  mk 
measurement  vectors:  Zk  =  {Cfc(*)}2=V  This  is  done  using 
models  for  the  likelihood  p(Zk\xk)  and  the  transition  prior 
p{xk\xk—i). 

2.1  Kinematics 

The  kinematic  state  xk  will  typically  contain  position  and 
velocity,  and  possibly  accelerations,  heading,  maneuver 
strengths  etc.  For  easy  of  presentation  we  assume  a  linear 
kinematic  transition  prior, 

xk  =  Fxk_ i  +  vk,  vk  ~  J\f(0,  Q).  (2) 

The  tracking  literature  provides  no  established  model  for  the 
development  of  the  mean  target  power  dk  through  time.  As 
the  fluctuations  of  signal  strength  from  scan  to  scan  can  be 
very  severe,  such  fluctuations  should  be  accounted  for  in  the 
measurement  model  and  not  in  the  process  model.  It  is  very 
difficult  to  estimate  a  time-varying  target  power  accurately 
for  fluctuating  targets  [14].  It  makes  more  sense  to  talk  about 
the  target  power  as  something  being  constant  over  several 
(say  >  20)  data  frames.  Therefore  we  treat  the  target  power 
state  as  a  constant  parameter  and  write  dk  =  dk-\  =  d. 

2.2  Measurement  model 

We  parameterize  the  measurement  vector  into  a  kinematic 
part  zk(i),  an  amplitude  part  ak(i)  and  a  background  part 
qk(i )  so  that 

Cfc  (*)  =  ( i),ak{i),ql  (i)]T.  (3) 

In  order  to  make  the  estimation  problem  well-posed  we  as¬ 
sume  that  at  most  one  measurement  C,k  ( i )  can  originate  from 
the  target.  All  the  other  measurements  are  thus  considered 
as  false  alarms.  We  also  assume  that  the  mapping  x  — >  2 
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is  linear,  for  example  through  application  of  the  coordinate 
transform  described  in  [4], 

zk  =  Hxk  +  wk,  wk  ~7V(0,  R).  (4) 

The  amplitudes  ak(i)  are  under  the  hypotheses  of  noise  or  a 
target  embedded  in  noise  assumed  to  be  distributed  accord¬ 
ing  to  likelihoods  p0{ak{i)  \  gk(i))  or  pr(afc(i)  |  dk,qk(i )) 
respectively.  These  will  be  thoroughly  discussed  in  Sections 
3.6  and  3.7. 

2.3  The  CLT  and  the  Rayleigh  model 

The  background  model  most  often  encountered  in  the  track¬ 
ing  literature  is  the  Rayleigh  distribution,  which  is  valid  un¬ 
der  the  Central  Limit  Theorem  (CLT).  The  backscattered 
signal  zk  in  resolution  cell  i  has  under  the  CLT  a  complex 
Gaussian  distribution.  The  envelope  or  amplitude  a  is  then 
Rayleigh  distributed, 

a  =  \z\  ~  Ra(a ;  77)  =  ^cxp  ^  ^  ^  (5) 


This  has  motivated  researchers  to  introduce  the  compound 
A' -model  [15],  in  which  radar  sea  clutter  is  modeled  by 
treating  as  a  correlated  “texture”  component,  while  a\r]  is 
an  uncorrelated  “speckle”  component.  We  let  77  refer  to  both 
the  deterministic  Rayleigh  parameter  in  (5)  and  to  the  tex¬ 
ture  random  variable  in  (7).  The  context  will  make  it  clear 
how  77  is  supposed  to  be  interpreted. 


2.5  Model  for  target  plus  clutter 

The  compound  interpretation  of  the  A' -distribution  allows  us 
to  evaluate  the  PDF  of  a  Swerling  I  target  in  A'-distributed 
clutter  as  follows.  Conditioned  on  the  texture  77,  the  back¬ 
ground  “noise”  w  has  a  complex  zero-mean  Gaussian  PDF, 


p(w\rj)  =  Afc(v)  —  N 


Re(tu) 

Im(tu) 


0  ' 

’  77  0  1 

1 

0 

5 

O 

_ 1 

The  complex  backscatter  signal  s  from  a  Swerling  I  target 
with  power  d  also  is  Gaussian, 


p(s|ti)  =  Afc(d). 


where  77  is  the  power  of  the  background. 


2.4  The  A -distribution 


The  CLT  is  valid  under  the  assumption  of  a  large  and  non- 
random  number  of  independent  scatterers  [13].  In  reality  the 
number  of  significant  scatterers  per  resolution  cell  may  not 
be  high  enough  for  the  CLT  to  hold.  One  often  observes  a 
more  heavy-tailed  background  than  expected  according  to 
the  CLT.  The  possibly  most  popular  model  for  such  non- 
Rayleigh  clutter  is  the  A'-distribution  [15].  A  AT -distributed 
amplitude  a  has  the  PDF 


KPDF  (a ;  77,  6)  4  - Kv_y  f  .  (6) 

Vb  r(i/)  \VbJ 

For  v  — »  00  the  K -distribution  turns  into  the  Rayleigh  dis¬ 
tribution,  while  v  <  1  indicates  a  very  heavy-tailed  back¬ 
ground.  By  r(-)  we  denote  the  Gamma  function  while 
Kv  ( • )  refers  to  the  modified  Bessel  function  of  the  second 
kind  [10]. 

The  A'-distribution  is  popular  for  two  reasons:  Physical 
plausibility  and  mathematical  convenience.  The  first  reason 
has  to  do  with  the  fact  that  reasoning  similar  to  the  CLT 
leads  to  the  A'-distribution  [2,15], 

In  contrast  to  so-called  stable  distributions,  the  A'- 
distribution  has  finite  moments. 


mn  =  E  [«”] 


_/-nr(i  +  f)r(z7+t) 

rM 


It  is  also  convenient  that  the  Af-distribution  can  be  viewed 
as  a  Rayleigh  distribution  Ra(a ;  77)  modulated  by  a  Gamma 
distribution  Ga(?7 ;  77, 6/2).  In  mathematical  terms. 


p(v ;  G  6)  =  Ga(r? ;  77, 6/2)  = 


Vv  1exp  (— ^y) 
(6/2)"r(77) 


(7) 


Therefore  the  sum  z  =  s+w  of  signal  and  noise  conditioned 
on  77  is  Gaussian  as  well,  and  the  corresponding  amplitude 
a  =  \z\  is  Rayleigh, 

p(z\d,  77)  =  A fc(d  +  77)  =>  p(a\d,  rj)  =  Ra(d  +  77). 


Under  the  compound  A'-model  the  texture  77  is  random,  so 
the  PDF  of  the  amplitude  must  be  found  by  marginalizing 
over  77, 


P(a\d,  A  6)  : 
a 

~  6^r(77) 


Ra(a;  d  +  if  dr] 


00  t?1'-1 


i]  +  d 


exp 


q2  ) 
2(77  +  d)  J 


dr]. 


(8) 


3  Methodology 

Most  tracking  systems  employ  a  modularized  architecture 
in  which  a  detector  extracts  measurements  to  be  fed  to  the 
data  association  and  state  estimation  methods.  This  does 
not  necessarily  have  to  be  the  case,  as  exemplified  by  TBD 
methods  which  do  not  carry  out  any  measurement  extraction 
at  all.  Although  our  work  is  relevant  for  such  methods  as 
well  [5],  we  are  in  this  paper  content  with  the  modularized 
architecture  as  represented  by  the  PDAF  and  the  PDAFAI 
methods. 


3.1  Extraction  for  Rayleigh  case 

For  any  cell  i  a  detection  is  declared  if  its  amplitude  exceeds 
a  certain  threshold  TlRk.  The  detection  threshold  TRk  is  de¬ 
termined  from  a  set  of  M  auxiliary  cells  (J'R  according  to 
the  closed  form  formulas 


% 


T\ 


Rk 


1 

2 M 


E  «>2 


%■ 


(9) 


p{a\rf  =Ra(o;  77), 

=>  p(a ;  77, 6)  =  KPDF  (a ;  77, 6) . 
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Here  fjk  is  the  Maximum  Likelihood  Estimator  (MLE)  of 
the  background  strength.  The  detector  given  by  (9)  and  (10) 
gives  an  actual  false  alarm  equal  to  the  design  false  alarm 
rate  l\A  as  long  as  the  background  is  Rayleigh.  It  is  there¬ 
fore  referred  to  as  a  CFAR  (Constant  False  Alarm  Rate)  de¬ 
tector  [9]. 

The  set  of  auxiliary  cells  is  chosen  to  lie  in  the  vicinity 
of  the  cell  under  test,  although  with  guardbands  to  prevent 
them  from  being  affected  by  signal  leakage. 

Each  extracted  measurement  vector  contains  the  location 
of  the  detected  cell,  as  well  as  the  value  of  the  cell  and  the 
clutter  estimate  from  its  auxiliary  cells.  As  discussed  in  Sec¬ 
tion  2.2,  range  and  bearing  can  be  converted  to  cartesian  co¬ 
ordinates  using  simple  formulas  [4],  The  full  measurement 
vector  of  (3)  is  under  the  assumption  of  a  Rayleigh  back¬ 
ground  parameterized  as 

C  (*)  =  [zT(i),a(i),fi(i)]T. 

Here  z(i)  =  [x(i),y(i)]T  is  the  centroid  of  resolution  cell 
number  i  in  the  cartesian  coordinate  system,  which  may  or 
may  not  be  related  to  the  target  through  a  model  such  as  (4). 

3.2  Extraction  for  K- case 

If  the  background  is  significantly  heavy-tailed,  the  detector 
given  by  (10)  becomes  increasingly  inadequate.  Instead  of 
attempting  to  estimate  the  local  background  power  //,  we 
will  estimate  the  /f -distribution  parameters  v  and  b  from  a 
somewhat  larger  set  of  N  auxiliary  cells,  and  then  set  the 
threshold  according  to  the  design  false  alarm  rate  for  a  K- 
distribution  with  these  parameters.  A  Method  of  Moments 
(MoM)  estimator  recommended  by  [2]  uses  the  first  and  the 
second  moments  of  the  data  given  by 

™i  =  ai and  m2  =  ^  K)2- 

j&G'k  i^GK 

This  estimator  attempts  to  solve  the  equation 


We  can  now  for  any  resolution  cell  i  determine  a  threshold 
TlK  according  to  the  criterium 

P(a  >  Tk)  «  1  -  Kv  ( 2-^f)  =  PFA.  (14) 

Vb  r(z/)  V  vo  / 

This  is  a  nonlinear  equation  in  Tk  which  must  be  solved 
numerically.  Instead  of  working  directly  with  the  cumulative 
distribution  function  as  in  (14),  we  recommend  using  the 
square  of  its  logarithm  for  numerical  stability. 

A  test  using  the  threshold  Tk  is  carried  out  in  order  to 
decide  whether  cell  i  should  be  stored  among  the  detected 
measurements.  When  a  /^-distributed  background  is  as¬ 
sumed  we  must  store  two  background  estimates  for  each 
measurement  vector.  In  addition  it  is  convenient  to  store 
the  detection  threshold  TK(i): 

<(*)  =[zT{i),a{i),i>{i),  b(i),TK{i)]T 
=[zT{i),a(i),q(i)]T. 

3.3  Probabilistic  Data  Association 

The  PDAF  [4]  is  perhaps  today’s  most  popular  tracking  al¬ 
gorithm  due  to  its  pragmatic  compromise  between  efficiency 
and  robustness.  It  is  a  suboptimal  algorithm  which  at  each 
time  step  collapses  the  target  state  posterior  PDF  into  a  sin¬ 
gle  Gaussian  which  then  is  propagated  to  the  next  time  step. 
Thus  the  past  information  about  the  target  at  time  step  k  can 
be  summarized  by 

p(xk\Z  )  ~  J\f  (xk  ai/e|fc— i,  Pk\k_  i).  (15) 

The  lumping  (cf.  Figure  1)  implied  by  (15)  is  done  by  ex¬ 
pressing  the  state  estimate  at  time  k  as  a  weighted  average 
of  the  prediction  xk\k_i  and  state  estimates  conditioned  on 
the  latest  measurements  zk(i).  This  leads  to  the  following 
Kalman  Filter-like  equations  for  prediction  and  measure¬ 
ment  update  of  the  state  estimate  xk\k  and  its  associated 
covariance  Pk\k: 


m2  4vV2{v) 

m\~  T(v+\Y 


Since  no  closed  form  solution  can  be  found  the  equation 
must  be  solved  numerically.  As  explained  in  [2]  a  second- 
order  approximation  can  be  found  as 


v  = 


1 

4 


(12) 


While  [2]  recommends  (12)  as  a  starting  point  for  the  nu¬ 
merical  search,  we  simply  use  (12)  as  our  estimator  of  v. 
Figure  2  shows  that  for  very  heavy-tailed  clutter  {v  «  0.1) 
the  estimator  (12)  has  a  significant  bias,  and  the  numerical 
search  is  necessary  in  order  to  obtain  an  unbiased  estimator. 
But  for  moderately  heavy-tailed  clutter  ( v  >  0.5)  the  nu¬ 
merical  search  appears  to  be  overkill.  Having  estimated  v , 
we  estimate  b  as 

b  =  m2/v.  (13) 


Xk\k—  1 


=  Fx 


k—  l\k—  1 5 


(16) 


Pk\k-1  ~  FPk-l\k-lFT  +  Q 

mk 

Xk\k  —  Xk\k  —  i  ~b  Kk  ^  '  pk ( <0 (i) 


Pk\k  =  Pk\k—1  —  (1  ~  Pk{fy)KkSk\k-lKk  +  Pk 


where 


Kk  =  Pk\k_1HT  Sk\k_1 

Sk\k-1  ~  P-  Pk\k—  lP-  “b  Pk 

vk(i)  =  zk  (?')  —  Hxk  |k_i  =  zk(i)  zk  |fc_i 

~ mk 

pk  =  Kk  ^2  l3k(i)vk(i)isk(i)T  -  ukuk 

_i—l 

mk 

i 'k  =  '22^k(i)^’k{i)- 

i= 1 


K! 
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Performance  of  MoM-estimators  for  u  =  0.1 


Performance  of  MoM-estimators  for  u  =  1 


Figure  2:  Histograms  illustrating  the  distributions  of  u  when  the  true  shape  parameter  is  0.1  or  1. 


The  index  i  ranges  over  all  mk  extracted  and  validated  mea¬ 
surements  inside  a  validation  gate  around  zk \k~\-  The  vali¬ 
dation  gate  is  defined  by  the  criterium 

(zfe(i)  -  zk \k-i)T  Sk1  ( zk(i )  -  zk\k-i)  <  g2, 


where  the  gate  threshold  g  typically  has  a  value  in  the  inter¬ 
val  between  3  and  10. 

Derivations  of  the  association  weights  /3k(i )  can  be  found 
in  [4] .  In  the  conventional  case  without  AI  they  are  written 


e(0 


Pk(i)  = 


_  )  b+^2jJ1e(j) 

_ b _ 

W-EjL'i  e(f) 


i  =  l,...,mk 
i  =  0 


(17) 


where 


e(i)  =  e-(zkii)-^k\k-i)TSk1(zk(i)-zk{k_1)/2 


The  value  of  b  depends  on  the  model  for  the  false  alarm  point 
process.  The  parametric  version  of  the  PDAF  assumes  that 
the  number  of  false  alarms  has  a  Poisson  distribution  with 
parameter  XV,  where  V  is  the  volume  of  the  validation  gate 
and  A  is  the  spatial  clutter  density.  It  can  then  be  shown 
that  [4] 


b  =  (2Tr/g)^XV(l  -  PdPg)/(c^Pd). 

For  the  more  robust  non-parametric  version  of  the  PDAF  we 
replace  XV  by  mk,  thereby  yielding 

b=  (2tt  /  g)^  mk(l  —  PdPg) /{Cu^Pd)- 

3.4  PDAF  with  Amplitude  Information 

In  [11]  it  was  suggested  to  improve  the  performance  of  the 
PDAF  algorithm  by  using  the  amplitude  of  the  measure¬ 
ments  as  part  of  the  tracking  algorithm.  The  amplitude  in¬ 
formation  only  affects  the  association  weights  and  not  the 


state  estimation  itself,  as  illustrated  in  Figure  1 .  The  modi¬ 
fied  association  weights  are  given  by 

e(i)  =  e-(zV‘)-^klk-i)TSk1(zk(i)-£klk_1)/2 

■  [pi{ak(i)\dk,qk(i))  /  PD\  /  [p0{ak(i)\qk(i))  /  PFA] 

A  g— (zfcW— 2fc|fc-i)  Sk  1(zfc(*)-Zfc|fc-i)/2 

Kak  (a)  I  dk  i  qk  (())  i 

where  the  amplitude  likelihood  ratio  l(ak(i)\dk,  qk{i))  de¬ 
pends  on  the  observed  amplitude  ak(i),  the  target  power  es¬ 
timate  dk  and  a  description  of  the  background  which  is  con¬ 
tained  in  qk{i)- 

The  likelihoods  p0(. . .)  and  pi(. . .)  are  the  PDF’s  of  the 
amplitude  under  the  two  hypotheses  of  only  clutter  and  of 
clutter  plus  a  target  with  power  dk.  Their  exact  expressions 
depend  on  how  the  background  is  modeled.  Notice  that 
Po(. . .)  and  pi(. . .)  are  divided  by  Pfa  and  Pd  respectively 
in  order  to  compensate  for  the  fact  that  ak{i)  has  exceeded 
the  detection  threshold. 


3.5  The  target  power 

Under  the  Rayleigh  assumption  an  MoM  estimator  dk  of  the 
mean  target  power  d  at  time  step  k  can  be  written 

4=max(A,A_Tt)  (,8, 

where 


Bk  = 


^  /  L  ^  mk_i- |_i 


2  L 


E 


\l  =  1  mk  /+1  2=1 

mk-i+i 


E  ak-l+l(i)flk-l  +  l(i) 


1  /  -  ^  i  ^  y 

Gk  =T  E ; -  E  Vk-i+i{i)Pk-i+i{i) 


L  \  z '  77Tju_/_i_i 

\Z=1  K  t  +  i  2=1 


B>k  —  £  £(1  -  /5fc— z+i (0)) . 

^  z= t 
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This  estimator,  which  is  developed  in  [7],  is  also  used 
to  estimate  the  target  power  when  the  background  is  K- 
distributed.  Then  is  found  as  half  the  second  mo¬ 

ment,  or  equivalently 


Vk-l+i(i) 


i>k-i+i(i)bk-i+i(i) 


2 


(19) 


The  time  lag  L  is  a  tuning  constant  whose  optimal  value 
hardly  can  be  determined  rigorously.  It  would  need  to  be 
rather  large  in  order  to  provide  an  accurate  estimate  of  cl.  On 
the  other  hand,  permanent  changes  in  the  target  amplitude 
may  in  practice  occur  quite  abruptly.  As  a  compromise  we 
have  found  L  =  20  to  be  a  reasonable  value.  When  the  track 
is  younger  than  L  time  steps  (i.e.  k  <  L)  we  necessarily  use 
a  shorter  lag. 

The  amplitude  estimate  plays  a  fundamental  role  in  the 
likelihood  ratio  introduced  in  Section  3.4.  It  is  also  useful 
in  order  to  determine  the  expected  probability  of  detection 
Pd-  For  the  case  of  a  Swerling  I  target  in  Rayleigh  clutter 
we  find  the  detection  probability  as 


threshold  (10).  In  the  context  of  TBD  this  has  been  shown 
to  have  a  cmcial  impact  on  performance  [5], 

The  estimator  77  given  by  (9)  is  a  random  variable  with  a 
corresponding  PDF  that  should  be  accounted  for  in  the  like¬ 
lihood  evaluation.  Its  PDF  p(fj\r])  depends  on  the  true  noise 
parameter  77.  To  utilize  the  information  carried  in  this  PDF 
requires  us  to  treat  77  as  random  as  well  and  assign  it  a  prior 
distribution  77(77).  By  marginalizing  the  joint  distribution  of 
a  and  77  given  77  we  obtain 


p{a\d,fi)  oc 


p(a\d,  v)p(y\v)p(v)dri- 


(22) 


The  background  estimate  77  given  by  (9)  is  under  the  as¬ 
sumption  of  independent  Rayleigh  samples  known  to  follow 
a  Gamma  distribution. 


77(77 1 77)  =  Ga 


?7M  xexp 

(&)"r(M) 


Since  we  do  not  have  any  more  prior  information  than  the 
value  of  77  itself  we  use  a  flat  and  thus  non-informative  prior 
in  the  Rayleigh  case: 


When  the  clutter  is  A' -distributed  we  find  the  detection  prob¬ 
ability  by  reasoning  analogous  to  (8)  as 


Pd  = 


(2 /by 

i» 


1/ —  1 

rj  exp 


T'k  \ 
2(77  +  d)J 


d?7  (20) 


77(77)  =  Uniform(77;  [0,  £]).  (23) 

Here  £  is  a  very  large  number.  The  flat  prior  allows  us  to  use 
77(77 1 77)  without  being  bothered  by  77(77).  The  likelihood  ratio 
becomes 


where  Yp  is  the  detection  threshold  as  obtained  from  (14). 
How  to  evaluate  (20)  efficiently  is  explained  in  the  Ap¬ 
pendix. 

3.6  Likelihoods  for  Rayleigh  case 

Under  the  Rayleigh-Swerling  I  assumption  used  in  [11]  the 
likelihoods  of  Section  3.4  become 


u  u  m  PpA  P(a\d^) 

l(a\d,r 7  =  — - 

PD  p(a|0,7?) 

=  Pfa  Jq  Ra  (a ;  77  +  d)  Ga  (77 ;  M,  jfc)  d 77 
pD  ft  Ra  (a ;  77)  Ga  (?)  ;  M,  ^ )  dr? 


(24) 


PoH,)  =  ?exp  (^) 

n(o|''-’,)  =  hT7i)exp(2STrf)) 

7/17  \  -Pfa  V  (  a2d  \ 

■—d-ap{^WTd))- 


(21) 


The  formula  (21)  is  optimal  given  that  the  true  values  of  d 
and  77  are  known.  Clearly  (21)  is  not  optimal  if  the  Rayleigh 
assumption  is  inadequate,  but  it  may  still  be  used  as  a  heuris¬ 
tic.  Several  papers  starting  with  [11]  have  demonstrated  the 
gain  from  (21)  under  ideal  circumstances.  The  performance 
of  (21)  under  less  ideal  circumstances  has  on  the  other  hand 
not  received  much  attention  in  the  tracking  literature. 

In  this  section  we  suggest  one  possible  strategy  of  cop¬ 
ing  with  the  inevitable  uncertainty  of  77  without  introducing 
the  A' -distribution.  This  strategy  is  to  account  for  the  uncer¬ 
tainty  of  77  in  a  way  analogous  to  what  is  done  in  the  CFAR 


Ka\d,v) 


Pfa  _  T(M) 

PD  [77M+  f]M 


( rj)M(v  +  d ) 


exp  - 


77  M 
77 


(25) 


q2  ) 

2  (v  +  d)J 


d/7. 


The  Appendix  explains  how  this  expression  can  be  evaluated 
in  an  efficient  way. 


3.7  Likelihood  for  K- case 

Significant  heavy-tailedness  is  under  the  compound  K- 
model  (7)  treated  by  assigning  a  Gamma  prior  distribution 
to  77.  The  uncertainty  imposed  on  77  this  way  is  likely  to 
be  more  severe  than  the  uncertainty  induced  by  the  CFAR 
background  estimate  (9).  Therefore  it  is  not  of  primary  im¬ 
portance  to  account  for  estimation  uncertainty  in  the  proper 
treatment  of  A'-distributcd  clutter.  Rather  we  simply  eval¬ 
uate  pi(a\d,u,b)  and  p0(a\v,b)  under  the  compound  K- 
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Pfa  —  0.005,  SNR  —  15dB 


Figure  3:  Various  likelihoods  for  SNR  =  15dB. 


model  (7)  as  done  in  (8).  This  yields 


Po{a) 

Pi(a) 
=>■  1(a) 


4a" 


VbV+1Y(v) 


AV-, [  — 


roo  ^-1 


6"r(i/) 

Pfa 


r]  +  d 


exP  T  ^ 


b  2(r)  +  d) 


d??, 


(aVb)1 


Pd  iK~-i(vs) 


f°°  ry"-1 


y  +  d 


V 


exp  ——  - 


b  2(77  +  d) 


dry.  (26) 


The  integral  in  (26)  must  be  solved  numerically,  as  outlined 
in  the  Appendix. 

Figure  3  shows  the  log  likelihood  ratios  In  Z(a| . . .)  corre¬ 
sponding  to  (21),  (25)  and  (26)  for  various  shape  parameters 
v.  Assuming  that  our  estimate  of  the  target  power  d  cor¬ 
responds  to  a  target  with  SNR  =  15dB,  the  curves  show 
the  log  likelihood  ratios  as  functions  of  the  actual  amplitude 
value.  Heavy-tailed  noise  is  seen  to  lead  to  a  more  conser¬ 
vative  likelihood.  For  really  spiky  clutter  (u  ss  0.1)  the  log 
likelihood  ratio  curve  corresponding  to  (26)  hardly  exceeds 
zero.  This  tells  us  that  tracking  a  15dB  target  in  such  noise 
most  likely  is  futile.  For  v  =  8  the  behavior  of  (26)  is  not 
very  different  from  (21),  which  means  that  the  Rayleigh  ap¬ 
proximation  may  be  considered  adequate.  The  behavior  of 
(25)  is  on  the  other  hand  always  qualitatively  different  from 
(21).  It  saturates  and  even  falls  back  for  very  large  a. 


4  Conclusion 

In  this  paper  the  classical  PDAF  method  has  been  tailored  to 
cope  with  heavy-tailed  clutter.  In  particular  we  have  devel¬ 
oped  conservative  likelihoods  which  evaluate  amplitude  in¬ 
formation  in  a  robust  way.  The  actual  performance  of  these 
developments  is  investigated  in  [6], 


Appendix:  Likelihood  evaluations 

The  tracking  methods  developed  in  this  paper  rely  on  eval¬ 
uations  of  the  integrals  in  (20),  (25)  and  (26).  Although  no 
closed-form  solutions  exist  to  any  of  these  integrals,  they 
can  be  evaluated  quite  fast  within  a  reasonable  accuracy. 
How  this  can  be  done  is  explained  in  this  appendix. 

According  to  [12],  the  integral  in  (26)  is  “readily  inte¬ 
grated  numerically”.  However,  the  integrand  has  several  un¬ 
desirable  properties  which  must  be  addressed.  On  the  one 
hand  it  may  tend  quickly  towards  00  as  77  — >  0,  while  on  the 
other  hand  it  may  tend  very  slowly  towards  0  as  ?y  — >  00. 
Its  overall  shape  is  difficult  to  characterize  in  simple  formu¬ 
las.  The  popular  strategy  of  importance  sampling  is  likely  to 
fail  since  the  obvious  factorization  may  lead  to  factors  with 
effectively  disjunct  supports. 

To  overcome  these  challenges  we  use  the  substitution 


7?  =  —  =>  dr/  = 
ir 


and  express  the  likelihood  ratio  as 


(27) 


l(a\d1  v ,  b) 


Pfa  _  (■ aVb 

Pv  2K-i  (Ji) 


m1-2"  ( 

it ^exp  r 


2 

u2b 


(28) 


u2a 2  \ 

1  +  u2d  J 


du. 


Denoting  the  two  integrands  g(rf)  and  g*(u ),  we  see  in  Fig¬ 
ure  4  that  the  behavior  of  g*  (u)  is  more  regular  than  the  be¬ 
havior  of  g(rj).  The  modified  integrand  g*(u)  is  integrated 
numerically  over  two  grids.  A  lower  grid  with  fixed  resolu¬ 
tion  captures  the  behavior  of  g*  ( u )  around  its  peak,  which  is 
given  by  a  cubic  equation.  An  upper  grid  with  exponentially 
decreasing  resolution  captures  its  asymptotic  behavior.  Ex¬ 
act  formulas  for  the  grids  are  given  in  [7],  but  will  be  left  out 
here  for  the  sake  of  brevity. 

Using  the  substitution  (27)  in  (20)  we  express  the  K- 
Swerling  I  detection  probability  as 


2(2/6)"  f°° 

r(z/)  J0  zt2"+1 


Tt2 


2(1  +  u2d) 


dit. 


This  modified  integrand  has  a  behavior  similar  to  g*  ( u ),  and 
the  same  technique  with  two  sampling  grids  can  be  used. 

The  substitution  (27)  also  makes  the  integral  in  (25)  more 
well-behaved: 


l(a\d}r)) 


Pfa  _  T(M) 
[f)M+  f]M 


,2M  — 1 


1  +  u2d 


exp 


r)M  + 


2(1  +  u2d) 


dzt. 


For  this  integral  the  grid  is  constructed  by  placing  100  sam¬ 
ple  points  uniformly  on  the  interval  u  £  [0,  yj\  / 77] . 
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Integrands  as  functions  of  rj  Integrands  as  functions  of  u 


Figure  4:  Example  integrands  of  the  A'-Swerling  I  likelihood  plotted  as  functions  of  the  texture  variable  rj  and  u  =  1  /  y/fj. 
The  amplitude  is  a  =  5  and  the  expected  SNR  is  15dB  for  all  three  cases.  The  scale  parameter  b  is  1  jv  in  order  to  keep  the 
second  moment  of  the  background  equal  to  one. 
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